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Abstract The primordial non-Gaussianity (PNG) in matter density perturbation is a very 
powerful probe of the physics of the very early Universe. The local PNG can induce 
a distinct scale-dependent bias on the large scale structure distribution of galaxies and 
quasars, which could be used for consnaining it. We study the detection limits on PNG 
from the surveys of the LAMOST telescope. The cases of the main galaxy survey, the 
luminous red galaxy (LRG) survey, and the quasar survey of different magnitude limits 
are considered. We find that the MAIN1 sample (i.e. the main galaxy survey with one 
magnitude deeper than the SDSS main galaxy survey, or r < 18.8) could only provide 
very weak constraint on PNG. For the MAIN2 sample (r < 19.8) and the LRG survey, 
the 2cr(95.5%) limit on the PNG parameter / NL are |/ NL | < 145 and |/nl| < 114 
respectively, comparable to the current limit from cosmic microwave background (CMB) 
data. The quasar survey could provide much more stringent constraint, and we find that the 
2<7 limit for |/nl| is between 50 and 103, depending on the magnitude limit of the survey. 
With Planck-like priors on cosmological parameters, the quasar survey with g < 21.65 
would improve the constraints to |/nl| < 43 (2a). We also discuss the possib ility of 
further tightening the constraint by using the relative bias method proposed by ISeljakl 
(120081) . 
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1 INTRODUCTION 

The seed of large scale structure (LSS) is usually thought to be produced during a period of in- 
flation in the early universe. Models of the standard scenario of single field slow roll inflation 
generally predict that these primordial perturbations are Gaussian distribute d dMaldacenal 120031 
lAcquaviva et al.L 120031 : ICreimneilil 120031: iBartolo. Komatsu. & Matarresel |2004|). However, significan t 



primordial non -Gaussianity (PNG) may be produced by multi-field inflation ( Linde & MukhanovHl997h 
and reheating (lDvali.Gruzinqv.& Zaldarriagai 120041) . if the inflaton potential has some sharp features 
( Chen. Easther. & Liml 20071), or i f the inflation is driven by field with non-canonical kinetic terms 



( Chen. Huang. Kachru. Shiu , 20071). PNG is also pr edicted by alternative models to inflation, such as 



the ekpyrotic scenario (ILehners & Steinhardtll2008l) . Thus, non-Gaussianity in primordial density per 



turbations could be a a very powerful discriminator for the p hysical models of the very early Universe 
(IBartolo. Komatsu. & Matarresel [2004t iKomatsu et all 120091) . 
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The PNG in density perturbation is usually parametrized as 

$(x) = 0(x) + / NL (0 2 (x)-<0 2 )), (1) 

where $ is Bardeen's gauge-invariant potential and <f> is a Gaussian random field. The parameter /nl, 
which denotes the effect of the PNG, is typically of order 10~ 2 in standard inflation models (e.g., 
iGangui et al.ll 1994 . Non-linear transformation from primordial field fluctuation to observables generi- 
cally gives rise to /nl ~ 0(1). With cp ~ 10~ 5 , this produces an extremely small non-Gaussianity in 
observables, such as the cosmic microwave background (CMB) temperature anisotropy or LSS bispec- 
trum. 

The CMB anisot r opy s o far offe r s th e cleanest observ at ional probe of PNG 
dWang & Kamionkowskii l2000t IVerde et all |2001|) . lYadav & Wandelt I d2008l) recently claimed 
detection of a positive /nl at 99.5% significance (27 < /nl < 147 at 2a) from the Wilkinson 
Microwave Anisotropy Probe (WMAP) 3-year data. Analysis of the WMAP 5-year data using 
the bis pectrum also favors a slightly positive /nl, —9 < /nl < HI at 95.5% confidence level 
(CL.t dKomatsu et all 120081) . while the Minkowski functional method yields -178 < /nl < 64 
(IKomatsu et al. L 120081) . In the future, the Planck satellite should significant ly improve the const r aint o n 



/nl, which may eventually push the uncertainty in /nl to be around 5 dKomatsu & Spergel 1200 lb . 
Alternatively, one may use the LSS observations to search for the PNG. However, even if the initial 
fluctuations are Gaussian, during the gravitational growth of the LSS, non-linear effect also generates 
non-Gaussianity, which may swamp the primordial one. Typically, the LSS constraints on PNG 
using the bispectrum technique are weak dVerde et all l2Q0lb 1 coccimarro. Sefusatti & Zaldarriagai 
l2004t ISefusatti & Komatsul l2007b . One consequence of PNG in structure formation is to significantly 
enhance the number of the rare peaks of fluctuation, boosting the fo rmation of massive objects 
dLucchin. Matarrese & Vittoriol 1 1 9861: iMatarrese. Verde & Jimenea 20001). One could o bserve this ef- 



fect from the abundance of galaxy clusters ( Benson. Reichardt. & Kamionkowskii |2002|) or dark matter 



halos harboring the first generation of stars dChen et all 12003 ). However, the required observational 
task is difficul t due to low - numb er statistics. 

Recently, iDalal et all d2008l) proposed a powerful and practical method to probe the local PNG 
with LSS observations. The PNG-caused enhancement on the formation of massive dark matter halos 
induces a distincti ve scale-dependent bias on the largest scales. Subsequent works confirmed and ex- 
tended this result dMatarrese & Verdel 120081: ISlosar et all 120081: lAfshordi & Tolleyl 120081: iMcDonaldl 



120081 : ICarbone. Verde. & Matarresell2008HKamionkowski. Verde. & Jimenez , 2008b . showing that use 



ful limits on PNG could be obtained from LSS surveys. 

In the present work, we investigate the sensitivity to PNG from the extragalactic surveys of the 
Large sky Area Multi-Object Spectroscopic Telescope (LAMOSTf]- The LAMOST is a recently built 
Chinese 4-meter Schmidt telescope with a field of view of 20 deg 2 . Equipped with 4000 optical fibers 
that are individ ually positioned by com puter during observation, it can take spectra of 4000 targets 
simultaneously dciiul tl 9981: Istonel 12008b . The LAMOST is currently undergoing engineering tests and 
calibrations, and it is expected to begin conducting surveys by 2010/2011. In a recent paper, we have 
forecasted the cons training power of the LAMOST extragalactic surveys on dark energy measurement 
dWang et all 12008b by considering different samples of galaxies and quasars. Here, we make similar 
forecast on the constraining power on PNG. 

2 METHOD 

On large scales, the PNG induces a distinct scale-dependent bias in the distribution of tracer objects, 
which can be measur ed from the correlation function or p ower spectrum of galaxies and quasars. 
Extending the work of IDalal et all d2008b . ISlosar et al ] d2008b found that this bias departures from the 
Gaussian case as 
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where b is the Eulerian bias and can be calculated as usual with the halo model, S c = 1.686, D(z) 
is the linear growth factor normalized to be the scale factor a in the matter domination epoch, O m o is 
the present matter energy density parameter, T(k) is the transfer function, Ho is the current Hubble 
parameter and c is the speed of light. 

The parameter p is introduced to account for the effect of merger bias. If the effect of merger on 
the formation of tracer objects is negligible, its value should be unity. However, the formation of some 
objects may be trigger by mergers. For example, quasars is believed to be accreting massive black holes 
in the center of galaxies, and the activity of quasar may be triggered by galaxy mergers, which send gas 
to the center of the galaxy, providing the materials needed for accretion. If the merger effect is dominant, 
p = 1.6. In general, the value of p is between 1 and 1.6. 

We consider three types of extragalactic surveys for LAMOST. (1) The main galaxy survey, which 
is a magnitude-limited general survey of all types of galaxies. By MAIN1 and MAIN2, we denote 
respectively the samples which are 1 and 2 magnitude deeper than the SDSS surveys, i.e. MAIN1 for 
r < 18.8, and MAIN2 for r < 19.8. (2) The LRG survey, a survey of color- selected luminous red 
galaxies (LRGs). Here for illustration we have used the MegaZ LRG sample dCollister et. all 120061) 
for calculation, though in reality the LAMOST LRG will be selected independently. (3) The quasar 
survey. Depending on the mag nitude limit, we defi ne three samples, QSOl (g < 20.5), QS02 (g < 21), 
and QS03 (g < 21.65). As in I Wang et all d2008l) . we assume that the surveyed area is 8000 deg 2 , for 
which we have the SDSS imaging data, even though the actual survey may cover a larger area if input 
catalogue for the south Galactic cap is available. For more details of these samples, including their 
number densities, redshift distributions, and the values of bias (estimated in the Gaussian models), see 
IWang et al.l (120081) . 

We model the galaxy power spectrum as 



P 9 (£0 = (6 9 + A6) 2 P lin (fc)i±^, (3) 

where b g is the bias of the galaxies calculated for the Gaussian case, Pu n is the linear matter power 
spectrum dEisenstein & Hul fl99&). k is the wavenumber (10~ 2 hMpc" 1 < k < 0.2 /iMpc -1 ), Q and 



A are free parameters which describe the non-linear effect. For fidu cial models, we fix Q = 4.6 for 
the m ain samples, Q = 30 for the LRG sample, and adopt A = 1.4 dCole et all 120051 : iTegmark et all 
120061) . For the quasar power spectrum, we adopt the simpler form 

P q (k) = (6 g + A6) 2 P lin (fc), (4) 

where b q is the bias of quasars. The value of bias for each of the sample has been given in lWang et al.l 
d2008l) . We assume that the PNG is small, so that the value of bias does not change much. 

The deviation of bias could be discovered with precision measurement of the power spectrum 
of the galaxy or quasar distribution. The mea surement error of the power spectrum is given by 
dFeldman.Kaiser.Peacocki 1 1 9941 iTegmarkl 1 1 997b 



fjp / 1 1 + nP _ j 1 

~P ~ V V s k^Ak nP ~ V V cS k* Ak ' (5) 



where V s is the survey volume and V e g is the effective volume, defined by 



off 



n(r)P 



n(r)P + 1 



2 

dr = 



nP 



n 2 



nP+ 1 



V s . (6) 



The second equality holds if th e comoving numbe r density n(r) is a constant. We divide the samples 
into different redshift bins as in IWang et al l d2008l) . then the total likelihood function is the product of 
the likelihood function for all the different bins. 

In order to estimate the error on the measurement of /nl, we use a Markov Chain Monte Carlo 
(MCMC) technique to explore the parameters space, then marginalize over (i.e., integrate out) the other 
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Fig. 1 Fiducial values of bias and the fits with the functional form b(z) = bo(l + z) 1 for 
the Mainl, LRG, and Q SQ3 samples. For each sample, the fiducial bias values (symbols) are 
from lWang etaDHooi, the thick curve are the best fit with the thin curves denotes the central 
68% distribution maginalized over the MCMC results. 



param eters to obtain a const raint on the PNG. For details of the MCMC technique, see, e.g., our earlier 
paper dGong & ChenL l2007h . The fiducial cosmology we adopt is the WMAP 5-year best fit ACDM 
model without PNG, with Q b0 = 0.05, Q m o = 0.24, fl A0 = 0.76, a s = 0.82, n s = 0.96, /nl = 0.0 
and ho = 0.7. The allowed ranges of parameters in the Monte Carlo are fibo £ (0, 0.1), Q m o £ (0, 1), 
(T 8 £ (0.6, 0.9), n s £ (0.9, 1.1), / NL £ (-400, 400), Q £ (0, 40) and h £ (0.5, 0.9). 

The value of bias for the sample is an important parameter which is partially degenerate with /nl- 
We model the value of bias at different redshifts bins in two possible forms: 

b(z) = b Q (i + zr, (7) 

or 

b(z) = a + c 2 /D (z) (8) 

where Do(z) is the growth factor normalized such that Dq(z = 0) = 1. This latter form is inspired by 
b(z) = 1 + (bp — l)/DQ(z), which can be derived in the linear regime for a passively evolved tracer 
population without any mergers by using S g = S m and 5 g {z) = b(z)S m {z) JF^ [l996Y Here b p is the 
present day bias factor of the tracer, S g and 5 m are the tracer and matter density contrast, respectively. 
We then try to use these two parametrization forms to fit the bias for our samples which were given in 
IWang et ail (120081) . Our test show that for the second form (Eqj8]l to yield good fit, we must have c\ < 0, 
which is physically unplausible given the original motivation for this form of bias evolution. Thus, we 
choose to use the first form (Eq. |7J to fit for biases at different redshifts. The parameter bo and 7 are 
then free parameters in the MCMC with the priors given by 60 £ (0, 10) and 7 £ (—5, 5). The predicted 
va lues of bias given b y the best fit parameter for several samples (Mainl, LRG, QS03) and those given 
bv lWang et aTTd2008l) " are shown in Fig.Q] We see that the functional form we adopt (EqO is adequate 



to describe the value of bias as a function of redshift for each sample. 

Note that to see how the LSS probes the PNG, here we do not consider joint constraint on /nl from 
CMB observations (e.g., bispectrum). We also do not use the currently available CMB angular power 
spectrum data in our MCMC unless otherwise stated (see next section about Planck prior). 
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3 RESULTS AND CONCLUSIONS 

The result of our analysis, i.e., the lcr(68.3%), 2er(95.5%), 3cr(99.7%) limits of /nl inferred from the 
power spectrum of the galaxy and quasar samples are given in Table Q] and Table [2] and also shown in 
Fig. [2] We see that for the MAIN1 sample, the 3cr(99.7%) limit on |/nl | reaches 379, so it provides only 
very weak cons traint on PNG. Th e MAIN2 and LRGs surveys have almost the same effective volume 
on large scales dWang et al l l2008h . so th ey provide compara ble constraint on /nl- Comparing with the 
current SDSS spectroscopic LRG survey dSlosar et al.ll2008b . the constraint on | /nl | could be improved 
by a factor o f about 1.5, to reach a 2er limit of about 1 14, which is comparable to current limits from 
WMAP data dKomatsu et aUl2008l) . 

Quasars are luminous and can be observed at high redshifts, so one could obtain large survey volume 
with qu asars. For the large scales considered here, the quasars provide the largest effective volume (c.f. 
Fig. 3 of lWang et al.ll2008t) . so we expect the LAMOST quasar surveys to provide a better constraint on 
/nl- To account for the possible effect of mergers, we calculated for the case of p = 1 (non-merger) 
and p = 1.6 (merger). We see that even for the QSOl sample (g < 20.5), the 95.5%(2er) limit on |/nl 
reaches 103 for the merger case and 82 for the non-merger case. For the QS02 sample (g < 21), the 
corresponding limits are 82(merger) and 63 (non-merger). For the QS03 sample (g < 21.6 5), the limits 
shrink to 69 (merger) and 50 (non-merger). Comparing with the SDSS quasar surveys dSlosar et al.1 
120081) . the LAMOST c ould improve the c onstraint of /nl by a factor of 1.3. We note that we assumed 
a lower QSO bias than lSlosar et al. I (l2008l) . so the limits obtained here are conservative. 

Table 1 The ler, 2a, and 3cr range of /nl constrained with the L AMOST Galax y Surveys. 
The last row (LRG-a) is from the relative bias method proposed bv lSef)akl(l2008h . 



Survey 


la 


2(7 


3(7 


Mainl 


±127 


±252 


±379 


Main2 


±71 


±145 


±206 


LRG 


±53 


±114 


±153 


LRG-a 


±13 


±28 


±37 



Table 2 The la, 2a, and 3a range of /nl constrained with the LAMOST Quasar Surveys. 



Survey 


Merger (p = 1.6) 


Non-merger (p = 1) 


1(7 


2a 


3(7 


1(7 


2(7 


3(7 


QSOl 


±52 


±103 


±155 


±40 


±82 


±122 


QS02 


±39 


±82 


±118 


±31 


±63 


±95 


QSQ3 


±33 


±69 


±101 


±25 


±50 


±75 



The posterior probability distribution function of /nl appears to be nearly Gaussian here, so the la, 
2a, and 3a errors are almost uniformly spaced. We have also investigated the degeneracy between /nl 
and ft m .o and n s . We find th at /nl is positively correlated with ft m .o, but it is almost independent with 
n s , confirming the results of ISlosar et al.l (120081) . 

By the time that the LAMOST surveys are completed, the Planck satellite should have been taken 
CMB for a few years. With these CMB data, the uncertainties on cosmological parameters should shrink 
further. Here we simulate this effect by adding priors on the other cosmological parameters based on 
the forecast of the Planck surveys. We set the standard deviation as fibo = 0.05 ± 4.0 x 10~ 4 , fl m o = 
0.24 ± 5.0 x 1Q- 3 , a H = 0.82 ± 6.7 x 10~ 3 , n s = 0.96 ± 4.0 x 10~ 3 and h = 0.7 ± 6.9 x 10" 3 
(IColombo et al.l 12008b . We present the results for the LRG survey and QS03 survey, and the impacts on 
other cases are similar. With the stronger priors, the constraints are improved slightly. With the QS03 
sample, we obtain |/nl| < 65 (merger case) and |/nl| < 43 (non-merger case) at 95.5% C.L., in 
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Table 3 The la, 2a, and 3a constraints on /nl from the LAMOST surveys with Planck 
priors. 



Survey 


la 


2(7 


3(7 


LRG 


±34 


±69 


±102 


QS03 merger 


±32 


±65 


±93 


QS03 non-merger 


±21 


±43 


±60 



Mainl 



Main2 - 
LRG 
LRG -a 
LRGa-r _ 
LRG+P 
QSOl 
QS02 
QS03 
QS03+P ■ 
QSOl-m- 
QS02-m- 
QS03-m- 
Q803-m+P- 



h±+ 
I 1 1 



I h 



I 1 h- 

I 1 h 

I 1 h 



H 1 



+-H 



H — h 



—\ 1 1 

H — I — I 



I I I ■ I I I 



II I ■ I I 



H h 



H 1 — 

-I — I — I 



H — I — I 



-200 



200 



Fig. 2 The la, 2a, and 3a range of /nl from LAMOST surveys. The samples are MAIN1 
(r < 18.8), MAIN2 (r < 19.8), LRG (the MegaZ sample), LRG-a (the constraint with the 
relative bias method using the LRG and the photometric samples), LRGa-r (the degraded 
relative bias method considering bias evolution), LRG+P (LRG with Planck priors), QSOl 
(QSO with g < 20.5, non-merger case), QS02 (QSO with g < 21, non-merger case), QS03 
(QSO with g < 21.65, non-merger case), QS03+P (QS03 sample, non-merger case, with 
Planck priors), QSOl-m (QSOl sample, merger case), QS02-m (QS02 sample, merger case), 
QS03-m (QS03 sample, merger case), QS03-m+P (QS03 sample, merger case, with Planck 
priors). The assumed survey area is 8000 deg 2 . 



comparison with | /nl | < 69 and | /nl | < 50 without the Planck priors. These results are given in Table 
3 and also shown in Figj2] 

Comparing directly the density fields of tracers with different biases in the same survey volume 
helps to eliminate the mea surement error due to cosmic variance an d break the degeneracy with other 
cosmological parameters (ISeljakL l2008t iMcDonald & SeliakL 120081). The relative bias a = 61/62 is 
proposed to be used as the estimator of PNG JSeljaki 120081: McDonald & Seliakl 120081) . Under the 
assumption that the two tracers have no stochasticity in them, and one tracer is unbiased (b=l), the 

constraint of /nl would be improved by a factor of (2 [(fiiPi)" 1 ± (r^i-^) 1 ] ) 1 2 ■ The feasibility 
of this method depends on the number density of the tracers. Here we consider the LRG as the high bias 
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tracer, and a sample of galaxies with photometric redshift as the low bias tracer. The av erage redshift 
of the LRGs is z f=a 0.4. Using the halo model method described in IWang et al.1 d2008l) . we find that 
for a sample of galaxies with good photometry up to r < 21.5, the average comoving density is n = 
8.11 x lCT 3 /i 3 Mpc" 3 , and the average bias is b = 1 .08. This sample of galaxies satisfies the requirement 
for the low bias tracer and has available SDSS photometry. With the relative bias method, we find that 
the limit on /nl obtainable with LRGs could be improved by a factor of about 4. 1 , with the la, 2a, and 
3a limits on |/nl| being 13, 28, and 37, respectively. 

In fact, since the evolution of the number density and clustering of the two galaxy samples could 
be different, the assumption of perfect correlation between the two tracer populations may not be valid, 
especially for the large volumes used for the non-Gaussinity test discussed here. This could introduce 
an additional error. The detailed study of this effect is beyond the scope of the present paper, we plan 
to investigate it in subsequent studies. For the problem discussed here, crude estimates show that in the 
worst case, it may cause a degration of the statistical power by as large as 50%. The LRGa-r in Fig. [2] 
shows the degraded precision. 

In summary, we forecast the constraints on the pri mordial non-Gauss ianity from the LAMOST sur- 
veys. Using the scale-dependent bias effect proposed by Dalai et al. ( 2008J), and based on the estimations 
of the galaxy and quasar surveys for the LAMOST dWang et all 20081) . we find that the galaxy surveys 
provide relatively weak constraint on the /nl parameter, but the QSO surveys could provide strong 
constraints, with a 2a limit on |/nl| somewhere between 50 to 103, depending on the sample size of 
the LAMOST QSO survey. After adding the Plank priors for other cosmological parameters, the 0SO3 
surve y could give a 2a limit of 43 for |/nl|- We also considered using the relative bias method dSeljakl 
120081) . Correlating the LAMOST LRG sample with a photometric sample of galaxies with r < 21.5, 
for which b w 1 at z = 0.4, we obtain a 2a limit of |/ NL | < 28. These results show that the LAMOST 
could make useful constraints on the primordial non-Gaussinity. 
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